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Exploiting the peculiar properties of proximity-induced superconductivity on the surface of a 
topological insulator, we propose a device which allows the creation of a Majorana fermion inside 
the core of a pinned Abrikosov vortex. The relevant Bogolyubov-de Gennes equations are studied 
analytically. We demonstrate that in this system the zero-energy Majorana fermion state is separated 
by a large energy gap, of the order of the zero-temperature superconducting gap A, from a band 
of single-particle non-topological excitations. In other words, the Majorana fermion remains robust 
against thermal fluctuations, as long as the temperature remains substantially lower than the critical 
superconducting temperature. Experimentally, the Majorana state may be detected by measuring 
the tunneling differential conductance at the center of the Abrikosov vortex. In such an experiment, 
the Majorana state manifests itself as a zero-bias anomaly separated by a gap, of the order of A, 
from the contributions of the nontopological excitations. 
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I. INTRODUCTION 

A Majorana fermion is an unconventional quantum 
state with non-Abelian statistics. Until recently, the con- 
densed matter community viewed it only as a mathemati- 
cal tool designed to help solving some specific many-body 
problems, arising, for example, in the areas of the two- 
channel Kondo modeUr- 3 - and the quantum magnetism^. 

However, the study of topological quantum 
computing^ initiated a search for experiments where 
this state can be directly observed and manipulated. 
Several proposals have been put forward. They rely 
of a diverse set of systems: liquid helium^ 7 -, topologi- 
cal insulators (TI)£, superconducting heterostructures^, 
Px+zPy-wave superconductor s 10 ' 11 , non-centrosymmetric 
superconductor s 12 ' 13 , proximity-induced superconduc- 
tivity on the surface of T I 14 i 15 ; Refs. [l6| studied 
non-abclian topological orders and Majorana fermions 
in s-wave superfluids of ultracold fermionic atoms and 
also spin-singlet superconductors with the spin-orbit 
interactions. 

In this paper we discuss a Majorana state localized 
at the core of an Abrikosov vortex residing in a two- 
dimensional (2D) superconductor. Clearly, not every su- 
perconductor has such a state inside its Abrikosov vor- 
tices: ordinary s-wave superconductors, for example, do 
not. Yet, in the theoretical literature several super- 
conducting systems are discussed where a vortex can 
trap a Majorana stat o 10 ' 11 ! 17 . However, these proposals 
have one serious drawback: in addition to the Majorana 
fermion, inside the normal core of the vortex, numer- 
ous non-topological Caroli-de Gennes-Matricon (CdGM) 
states are localized as well^. These states are separated 
from the zero energy by a minigap S whose size can be 
estimated as: 
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where A is the superconducting gap and is the Fermi 
energy. In order for a device to be a building block of 
a topological quantum computer it is necessary to freeze 
out all non-topological degrees of freedom, that is, the 
operational temperature should be much smaller than S. 
Since ej? > A, the minigap is expected to be extremely 
low (of the order of 10 -2 — 10 -3 K for usual s-wave su- 
perconductors). This means that such proposals have 
very dim prospects, unless a way of increasing 5 is found 
(however, see Ref. Il9f ). 

A possible way to overcome this shortcoming is de- 
scribed in Ref. [lH (see also Ref. [l6h . In this reference 
the idea of a "robust" Majorana fermions is put forward: 
the Majorana state is robust if the eigenenergy of the 
lowest non-topological excitation is of the order of A, 
that is d ~ A. It is found numerically that, if conditions 
are right, the robust Majorana fermion exists in a vortex 
residing in a proximity-induced superconductor on the 
surface of a TI. This result implies that the Majorana 
fermion in such a system can be created and manipu- 
lated at experimentally achievable temperatures. 

A device suitable for this task is presented in Ref. [ID 
as well. It relies on two coupled tri-j unction devices in 
which a Josephson vortex is inserted. A tri-junction is 
a meeting point of three Josephson junctions separating 
three superconducting islands placed on the surface of 
the TI. Altogether, the system consists of four super- 
conducting islands and four superconducting loops with 
magnetic fluxes to control the superconducting phases 
on the islands. The Majorana fermion is bound to the 
Josephson vortex. Varying the relative phases with the 
help of the fluxes one can move the fermion from one 
tri-junction to another. 

Here we discuss a much simpler system in which the 
robust Majorana fermion may exist, as shown in Fig. [TJ 
It is related to the proposal of Ref. UH: the most basic 
component is the vortex inserted into the superconduc- 
tor induced on the TI surface by the proximity effect. 
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It is demonstrated below that such a vortex can host a 
robust Majorana state whose presence can be detected 
with the help of local tunneling experiments. Thus, our 
proposed setup can provide a proof-of-principle that a ro- 
bust Majorana fermion is indeed possible and robust, as 
claimed. However, the simplicity comes at a price: unlike 
the device of Ref. Il5l . our Majorana fermion is pinned in 
space. 

In addition to that, our results are as follows. We 
investigate the Bogolyubov-de Gennes (BdG) equations 
which describe an Abrikosov vortex in our system and 
analytically demonstrate that, indeed, the gap, separat- 
ing the Majorana state and the lowest non-topological 
excitation is of the order of A, in agreement with nu- 
merical results^. Further, we provide simple arguments 
explaining why the robustness of the Majorana fermion 
exists in this system: it is a consequence of the vanishing 
density of states in the TI. 

Our paper is organized as follows. In Section II we re- 
view the derivation of the BdG equations on the surface 
of the TI for the uniform case and derive BdG equations 
in the presence of an Abrikosov vortex. In Section III we 
analytically obtain the zero-energy solution (Majorana 
state) of these equations. In Section IV, excited states of 
the model and the robustness of the Majorana fermion 
are analyzed. In Section V we discuss the results ob- 
tained here. More technical points are relegated to two 
appendices. 



II. BOGOLYUBOV-DE GENNES EQUATIONS 
A. General formalism 

The system under investigation is schematically pre- 
sented in Fig.[TJ It consists of a TI sample on which a slab 
of s-wave superconducting (SC) material is placed. On 
the surface of the TI a 2D band of electron states exists. 
This band is described by the massless Dirac equation 
(Weyl-Dirac equation). The proximity to the supercon- 
ductor induces a finite gap Ati in the Dirac band. 

When such a system is placed into a transverse mag- 
netic field of sufficient magnitude an Abrikosov vortex 
enters it. This vortex is accompanied by a "pancake" 
vortex inside the 2D Dirac band of TI. Such a pancake 
vortex can host a Majorana stated. Since the core of 
the Abrikosov vortex contains a large amount of CdGM 
states separated by small energy gap 5, the Majorana 
state is not robust. To counteract this disadvantage a 
cylindrical channel of radius R > £ is carved in the su- 
perconductor (see Fig.rfJ). The cavity removes the CdGM 
states from the vortex core. It also acts as a pinning cen- 
ter for the vortex. 

When the CdGM states are absent, the remaining low- 
lying states may be present only inside the "pancake" 
vortex core in the TI. To find them we derive the effective 
BdG equations for the TI degrees of freedom. To this end 
we consider the Hamiltoniani^ describing the proximity 



effect at the TI-SC interface 



H = H' 



TI 
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where Hti, H^c ar e the Hamiltonians for the TI sur- 
face and the BCS s-wave superconductor, and T (T') 
accounts for the tunneling from the TI surface to the SC 
(from the SC to the TI surface) . The excitation spectrum 
of the model is described by the equation 



(3) 
(4) 



where i?Ti,SC an d T are written as a 4 x 4 matrix in the 
Nambu basis^ 



H T1 = [iva ■ V r - U(r)] t, 
Hsc 
f 



2m 
r z T(r -!•'). 



f )t z + A'(R)r x 



(5) 

A"(R)r„(6) 
(7) 



In these equations h = 1; R = (x, y, z) is the 3D coordi- 
nate inside the SC; r = (x, y) is the 2D coordinate on the 
surface of the TI; a and r are the Pauli matrices acting 
in the spin and charge spaces, respectively. The param- 
eter v is the effective electron velocity at the TI surface; 
ep is the Fermi energy in the SC. The Fermi level U(r) 
in the topological insulator may be inhomogeneous: it 
depends on the external potential and the tunneling op- 
erator (see Appendix [A|. The tunneling kernel T(r — r') 
is independent of spin and charge indices. 

The wave functions "Jti.sc ar( 3 the 4-componcnt 
spinors: 



TLSC = [ut> u l' V l>~ v tY 



(8) 



The spinor ^ti = ^Ti(r) corresponds to the surface state 
and depends on x and y only. The spinor ^>sc = v I / sc(R) 
describes electrons in the superconductor bulk. It van- 
ishes for z < 0. 

The complex order parameter in the SC is A = A' + 
iA" , where both A' and A" are real. The superconductor 
is characterized by the correlation length 



VF 



(9) 



Starting from Eq. ((5|) and Eq. (j6)), the effective BdG equa- 
tions for the TI states can be derived^ 



H T1 - ft (H SC - w) _1 fj *ti = 



'TI- 



(10) 



Here the expression (i?sc — UJ )~ 1 is the single-electron 
Green's function for the superconductor. For details, see 
Appendix [A] 
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FIG. 1: (Color online) Schematic side view of the system 
studied in this work. On the surface of the topological insu- 
lator (TI) there exists a 2D band of surface states which may 
be described by massless Dirac equation. To induce supercon- 
ductivity in this 2D metal, a slab of superconducting material 
(SC) is placed on top of the topological insulator's surface. A 
cylindrical cavity of radius R is carved in the superconductor. 
It serves as a pinning site for the Abrikosov vortex, which is 
introduced in order to create the Majorana state. The main 
purpose of the cavity is to remove the unwanted single-particle 
excitations in the normal core of the vortex. Therefore, a 
robust Majorana fermion can be localized inside the vortex 
core created on the surface of the topological insulator by the 
proximity effect. 



where the effective Hamiltonian and its parameters are 
defined as 

H cff = [i«(w)((T • V r ) - £/■(«)] t z (12) 

+A'(u))t x + A"(w)r y , 
Uv /|A| 2 - 



V|A| 2 -w 2 + A' 



(13) 



^)= {u+ / ^ vl f7\ (14) 

y/\A\ 2 -LJ 2 + \ 

AH = A'H + iA» = /|A|2 AA 2 ■ (15) 

V |Ap - or + A 

The quantity A has the dimension of energy. It char- 
acterizes the transparency of the barrier separating the 
SC and TI, and can be measured in a tunneling experi- 
ment performed at T > T c . For our purposes we need a 
sufficiently thick insulating layer between the supercon- 
ductor and the TI to guarantee the low transparency of 
the barrier (A <C £f)^ 

The contact with the superconductor shifts the bare 
Fermi energy U by the amount SU = 0(A). The details 
of the derivation can be found in Appendix [A~] 

For uniform A and A we can choose A to be real. 
Eq. (fTTj) has no solutions with |w| < Ati, where the 
proximity-induced gap Ati < A satisfies the equation 
A(Ati) = Ati, or, cquivalently, 



A _ A TI I A + A T i 
A ~ 



A 



A- A 



TI 



(16) 



The gap in the TI is a monotonous function of A: Ati ~ A 
at small A <C A and approaches A from below at large 
A > A. 



B. Uniform system 



The symbolic Eq. (fTU)) is very general: it is valid for 
arbitrary A(R) and T, as long as the Green's function 
remains well-defined. Of course, finding the Green's func- 
tion for a non-uniform A may be practically challenging. 
Fortunately, in the regime of interest the required Green's 
function can be constructed from the knowledge of the 
Green's function for a homogeneous system. Thus, as a 
first step, let us study the situation when A and T are 
uniform over the whole interface, the cavity is absent, 
and \uj\ < \A\. The resultant BdG equation now reads 



C. The system with the pinned vortex 

The system schematically drawn in Fig. [TJ however, is 
not uniform. Due to the cavity and the vortex, both A 
and A acquire some coordinate dependence. Obviously, 
A(r) vanishes for r < R, where r = \/ x 2 + y 2 is distance 
to the axis of the cavity. In addition, when the vortex is 
introduced^, A(R) can be written as 

A(R)HA(r)|e ie , (17) 

where 9 is the polar angle in the (x, y)-plane, and |A(r)| 
is an increasing function of r, which approaches the bulk 
value |A| when r oo. It is finite at r = R + 0. In 
such a case, strictly speaking, one has to re-calculate the 
superconducting Green's function for a spatially varying 
A(r, 9). This might be particularly difficult for r < £, 
where the phase 9 varies quickly on the distances of the 
order of £. 

However, one can avoid the latter complication if 



(H ee - w) #ti = 0, 



(11) 



(18) 
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In this limit our formalism can be easily adopted to ac- 
count for the vortex presence. Ignoring the detailed be- 
havior of |A(r)| when r R, we assume that 



A. Majorana state for large R 



\A{r)\ = \A\0{r-R), 



(19) 



where i9(r) is the Heaviside step-function. 

When Eq. ([15} holds true, the order parameter phase 
8 varies slowly on distances of the order of £. Thus, it is 
permissible to insert the non- uniform A(r, 8), Eq. ([T7]). 
directly into Eqs. (fT2)) - (|T5|) . Since |A| is r-dependent, 
therefore, w, |A|, and U are non-uniform. Further, we 
assume that our treatment remains valid, at least quali- 
tatively, in the case R > £. 

In our formalism the TI area beneath the cavity (r < 
R) is non-superconducting. It may be viewed as the nor- 
mal core of "the pancake" vortex. Outside the core (for 
r > R), the absolute value of the order parameter equals 
to its equilibrium value |A|. This approximation is very 
natural in the case R > £. 

To calculate the eigenenergies of the Hamiltonian 
Eq. (|T2"j) with a vortex, it is standard to exploit the cylin- 
drical symmetry of the problem to separate the variables 
(for details see Appendix [Bj . If we define the spinor $ 
as 



*ti = exp [i8{t z + a z )/2 + ifx9]^(r), 
then, its four components satisfy 
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tf + W£-{L, + U)tf 



+ Z)ti + \A\tf-(w-U)fg 



|A|/ 2 " 



(20) 
(21) 



0, 
0, 

0,(22) 
0. 



The system of equations (|22|) can be solved exactly 
when ui = 0, fj, = 0, and SU(r) = 0. The latter re- 
quirement is satisfied only in the presence of the external 
gate potential which compensates for the Fermi level shift 
induced by the coupling to the superconductor (see Ap- 
pendix In subsection MI Bl it is demonstrated that a 
non-zero 5U (r) does not destroy the w = solution. 

To find the desired solution we define new functions 



x% = if? /i 1 , 



(23) 



For these functions the system of Eqs. (|22|) splits into two 
systems of two equations each: 



+ \A\X2 



dr 



dr 
\A\Y? 



= iUY? , 



-y" = iux. 



and 



v*f - |A|A? 
dr 



dY 2 ° 
i — — 
dr 



|A|F 2 + -y 2 ° 

r 



2 ■ 



iUY 2 ° , 



iUXl 



(24) 
(25) 

(26) 
(27) 



An elementary analysis reveals that the system of 
Eqs. (|21)1) - (|27D has no non-zero solution decaying at r — > 
co. Thus, X® — Y 2 a = 0. The non-zero solution of the 
system of Eqs. (|24l) - (f2"5|) can be written explicitly. Keep- 
ing in mind that U(r)/v(r) = U/v is a constant indepen- 
dent of r [see Eqs. (fT3"]l and (jl~4")) ]. one derives: 



Y? 



= 2C 



iJo(Ur/v) 
Ji(Ur/v) 



exp 



r dr'\K(r')\ 
v(r') 



(28) 

where C is a normalizing coefficient, and Jo, Ji are Bessel 
functions. Thus, the exact solution for uj = is 



This is the most general system of equations describing 
the sub-gap states near the cavity. Below we study the 
spectral properties of this system. 



III. ZERO-ENERGY MAJORANA FERMION 
SOLUTION 

In this section we will demonstrate that the Hamilto- 
nian H has a zero-energy Majorana fermion solution for 
arbitrary R. First, we will demonstrate this for the case 
of large R, when one can map H on H e g. Afterward this 
result will be generalized for any R > 0. This section 
contains some known results (see, e.g, Refs. [20ll22Tj24T ) . 
but these are included here to make the derivation more 
complete and self-contained. 



* M (r) = C 



( exp(i6)J (Ur/v) 
-Ui{Ur/v) 
MUr/v) 
\iexp(-i0)J o (Ur/v)y 



exp 



r dr'\A(r')\ 
v(r') 



(29) 

If the ratio \A(r)\/v(r) = X/v is a constant independent 
of r [sec Eqs. and (|15p] then the spinor decays 
for distances larger than v/X, which may be viewed as a 
characteristic localization length of the zero-energy state. 

To prove that the eigenfunction given by Eq. (|29[) cor- 
responds to the Majorana state, consider the following 
fermion operator 



<(r)^(r)+ U f(r)^( 



M / 



(30) 



+ Vt M (r)^ t (r)+z,f(r)^(r) 



Mi 
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where ipl (r) is the creation operator for an electron with 
spin a located at point r. The functions u^(r), v^(r) 
are components of the spinor ^ M (r). The operator 
creates a fermion in the state corresponding to $M(r). 
It is easy to demonstrate, by direct calculation, that 
= i& M - Therefore, V^m corresponds to the Majorana 
fermion. 



B. Majorana state for arbitrary R and U 

We demonstrated above that for large R our system 
can be mapped to 2D Dirac electrons. In the latter 
model, when the vortex is present, the zero-energy so- 
lution is found 2 ^. 

Unfortunately, if R is small, this mapping is inappli- 
cable. What happens to the Majorana state when the 
cavity is small? Below we will prove that our Hamilto- 
nian has a zero-energy eigenstate for any R > 0. 

We start our reasoning with the observation that H 
satisfies the following charge-conjugation relation: 



H 



-TyUyH -I 



yUy. 



(31) 



A. Analytical calculations 

To find the eigenfunctions for these states, it is neces- 
sary to solve the system (|22|) for generic w and \i. The 
solution can be simplified significantly for U = 0. We will 
now investigate this case. The non-zero U(r) may be ac- 
counted with the help of perturbation theory, at least for 
\U\ < A. According to Eqs. (26) and (27) of Ref.GjJ the 
case \U\ > A is not favorable for the robustness of the 
Majorana state and will not be studied here. 

When [7 = and r < R, the system (|22|) decouples 
into two sets of equations: 



d_ 

dr 



IV 
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d fx 
dr r 



o, 



0, 
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iv | 






K dr 
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' d 






dr 
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ft- "ft 



0. 



(32) 
(33) 

(34) 
(35) 



The solution is 



Thus, for every eigenstate \& of H with a non-zero 
eigenenergy u> 7^ 0, an eigenstate T y a y ^f* with eigenen- 
crgy —lu is present. A spinor with positive eigenenergy 
corresponds to the creation of a quasiparticle, while the 
charge-conjugated spinor corresponds to the destruction 
of this quasiparticle. 

Further, it is demonstrated here that for large R the 
Hamiltonian has the zero-energy solution This 
eigenstate is special for it remains unchanged after a 
charge-conjugation transformation. This means that the 
number of eigenenergies lying inside of the even-energy 
interval (|A|,— |A|) is odd (i.e., all the non-zero eigen- 
states are paired, while the Majorana state is unpaired). 
If we start decreasing R this property endures: due to 
symmetry [see Eq. (fBTj) ] the eigenstates can enter or leave 
our energy interval only in pairs. Thus, for any R the 
Hamiltonian H has an unpaired u> = eigenstate invari- 
ant under charge conjugation. However, it is necessary 
to remember that, when R < £, the Majorana state can- 
not be robust due to the CdGM states in the Abrikosov 
vortex core. 

Instead of R, one can vary SU. The above reasoning 
can be modified to prove that the deviation of SU(r) from 
the SU (r) = value does not destroy the Majorana state. 



IV. EXCITED STATES IN THE VORTEX CORE 

In addition to the cj = state, it is possible to have 
< \uj\ < IAtiI states localized at the vortex core. 



f^iB^-xiur/v), f£ = BMur/v), (36) 

where A M and are constants. 
For r > R the equations are: 



djq 

dr 



Al 



-xk 



v „_iv dX? |A 
1 ~ w dr 



d 2 Af 1 dX? 
1 + 



1 



dr 2 
dr 2 



dr 

ldX% 
r dr 



.KM] 2 
1 



|A| 

vr 



(37) 



[£( w )] 2 vr r 2 



where the functions A{' 2 and Yj_ 2 arc defined by 
Eqs. (f23|) . and the length £(w) is given by the formula: 

v{w) 



|A(w)| 2 - LO 2 



(38) 



The above equations may be solved 2 ^ in terms of the 
Whittaker functions W at p(z): 



A 1.2 - 



Q1 , 2 , M (2r/f(«)) 



<*1,2 = T- 



IAI 



Al 2 



(39) 
(40) 
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FIG. 2: (Color online) The energy gap 8 between the ground 
state and excited states as a function of the cavity radius R, 
for different transparencies A of the barrier. The solid curves 
1, 2, and 3 correspond to the states with n = 0, (i = 1, 
and the dashed curves a, b, c correspond to the states with 
n = 1, fj, = 0. As we go from top to bottom, the transparency 
A decreases: curves 1 and a (black) are drawn for A/ A = 2; 
curves 2 and b (red) are drawn for A/A = 0.5; curves 3 and c 
(blue) are drawn for A/ A = 0.1. 



Using Eq. (|16p one can show that for subgap states (|o;| < 
Ati) the expression under the square root in Eqs. (|38|) 
and (|4"0"j) is positive, therefore, and ax, 2 are real. 

Equation (|39[) implies that £(w) is the energy-dependent 
localization length for the subgap states. 

Matching the solutions at r = R we derive the equation 
for the subgap eigenenergies: 



B. Majorana fermion robustness 

The excited quantum states of our system can be clas- 
sified using two quantum numbers: the radial number n 
and the orbital number \i. The Majorana state obtained 
above corresponds to the n = 0, fx = state. The results 
of the numerical solution of Eq. (|41|1 for different values 
of (n, (j,) are shown in Fig. [2j Here the reduced gaps 
8/ A between the ground state and excited states with 
(n = 0, fj, = 1) and (n = 1, [i = 0) are plotted for differ- 
ent transparencies A as a function of the cavity radius R. 
The solution of Eq. (|4"Tj) confirms an intuitively transpar- 
ent conclusion: that the first excited state is either the 
(n = 0, /i = 1) state, or the (n = 1, fi = 0) state. Note, 
that, to calculate <5 in the case (n = 1, fj, = 0), we use the 
asymptotics of the Whittaker functions W Ql 2,0(2) valid 
at \z\ 3> 1. This negatively affects the accuracy of our 
calculations at small R. We believe, however, that this 
substitution does not distort the qualitative features of 
the solution. 

As it is seen from Fig. the state (n = 0, /i = 1) 
lies lower than the state (n = 1, fj, = 0). Thus, the 
gap between the ground state and the excited state 
(n = 0, /i = 1) characterizes the robustness of the Ma- 
jorana fermion in our system. The gap value increases 
when A increases. If A > A, the energy gap S(R) prac- 
tically saturates and further growth of the barrier trans- 
parency does not significantly improve the robustness. 
Experimentally, the regime A ~ A corresponds to a bar- 
rier with low transparency^, since A is much smaller than 
the Fermi energy. 

At a given A, the curve 5(R) is a decreasing function 
of R, approaching a maximum value 

5 max = A(0) = at (42) 

The main conclusion that follows from the results 
shown in Fig. [2] is that the energy gap between the Ma- 
jorana state and the first excited state may be of the 
order of A (8/ A = 0.5-0.6 or even higher) if A/A > 1-3 
and = 2-3. Therefore, a suitable choice of R and A 
allows one to realize the robust Majorana state. 
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(41) 



In this equation all functions W a ^(z) and W' a p(z) must 

be evaluated at z = 2R/£(w), and all Bessel functions 
must be evaluated at ujR/v. We analyzed Eq. (|4ip nu- 
merically. 



C. Physical and intuitive explanation of the 
Majorana state robustness 

Above we demonstrated that in our system a robust 
Majorana state exists. Our results agree with previ- 
ous numerical calculations^. However, it is desirable to 
have a simple non-technical physical argument explaining 
these results. To this end, now consider the "pancake" 
core. It can be approximately described as a circle of 
radius R, where the superconducting gap is zero. Let 
us now evaluate the number 7V£j re of the single-electron 
subgap states: 

/•A TI a2 c>2 

N^r-nR 2 VTi(e)de (43) 

Jo v 
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Here iAn(e) ~ e/v 2 is the density of states for the TI. The 
average energy interval between these states is 



>T1 



Ati/N. 



TI 



^ 2 /(i? 2 A TI ). 



(44) 



For R ~ £ and A T i ~ A, one has N%f° ~ 1 and S T i ~ A. 
The last estimate is equivalent to the statement of the 
Majorana fermion robustness. 

If, instead of TI, we now consider a 2D superconductor 
with parabolic dispersion, we then obtain: 



N c 



7ri? 2 Az/ m (e F ) 



v a) e 



(45) 



Here v m (e) is the density of states for a 2D metal [v m (s) — 
4m]. When R ~ £, the number A^™" is much greater 
than unity. Further, the energy difference between the 
CdGM energy levels is of the order of (A 2 /e F )(^ 2 /i? 2 ). 
For R ~ £ we recover Eq. ([1]). 

The estimate of the previous paragraphs demonstrates 
that the source of robustness of the Majorana fermion in 
TI is the vanishing density of states at the Dirac point. 
If we were to apply a bias shifting the Fermi level away 
from the Dirac point, at a bias value U substantially ex- 
ceeding A, we recover Eq. ([T]). Thus, U = is a favorable 
condition for the Majorana state robustness. 



V. DISCUSSION AND CONCLUSION 

We demonstrated that in the proposed system a ro- 
bust Majorana state may exist. A possible way to detect 
its presence is to perform the following tunneling experi- 
ment: insert an STM tip into the cavity and measure the 
electron current flowing through the TI surface. In such 
a setup the Majorana state manifests itself as a zero-bias 
anomaly: a peak of the differential conductance located 
at eV = 0. Unlike the zero-bias anomaly observed for an 
Abrikosov vortex in a BCS superconducto r 26 ^ 7 , in the 
case of the robust Majorana state, the zero-bias anomaly 
is accompanied by several discrete subgap peaks which 
correspond to the excited states bound in the vortex core. 

We would like to compare our proposal with the sys- 
tem discussed in Ref . HH. It is clear that in our system the 
Majorana fermion is pinned to the cavity. The advantage 
of the setup in Ref. [l5| is that it allows one to shift the 
Majorana fermion along a straight channel over the dis- 
tance L ~ £. However, the robustness of the Majorana 
state decreases when one tries to increase the distance L 
above £: the gap separating the Majorana fermion and 
the lowest excited state can be estimated as S ~ v/L for 
L > £. Moreover, the complexity of their— system is an 
additional limitation. 

To conclude, we propose a SC-TI setup in which a 
robust Majorana state may be realized inside the core 
of the pinned vortex. The robustness was justified with 
the help of both analytical and physically intuitive ar- 
guments. An experimental detection of this state is also 
discussed. 
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Appendix A: Derivation of the effective Hamiltonian 
for proximity-induced superconductivity 

In this Appendix we briefly present a derivation of the 
effective BdG equation at the TI surface. Here we refine 
the reasoning of Ref. [l5ll28l 

The system in question consists of a superconducting 
slab, described by the BCS Hamiltonian Hsc, and a topo- 
logical insulator with Hamiltonian Hti- These two are 
separated by a fiat interface. Tunneling across this inter- 
face is described by the tunneling Hamiltonian: T + ■ 

As the TI has the bulk gap, its low-energy states are lo- 
cated at the surface. At low energy these 2D states can be 
approximately described by the Weyl-Dirac Hamiltonian 
whose apex is located at the M point at the boundary of 
the Brillouin zone of the TI. 

Due to conservation of the spin and the quasi- 
momentum parallel to the flat interface, the tunneling 
Hamiltonian couples a single-electron state 



k<r) insidG 



XM+k,k x ,a) inside the 



the TI with a normal-metal state 

superconductor. Here a is the spin index, k is the mo- 
mentum's components parallel to the barrier measured 
from M, and k z is the absolute value of the trans- 
verse momentum (since scattering at the interface cou- 
ples SC states with k z and — k z , the state XM+k k, a) * s 
a boundary-condition-compatible linear combination of 
both k z and — k z states). The superscript 'm' stands for 
'metal'. Let us denote the corresponding matrix element 
by T(k, k z ). It is assumed to be spin independent. 

Invariance of the tunneling Hamiltonian with respect 
to the spatial inversion r — > — r implies that the same ma- 
trix element couples the states j^ 1 :^) and X'^M-k fc aj- 
Indeed, the quasi-momentum M — k of the TI state is 
equal to the momentum — M — k of the state inside the 
superconductor, modulo the TPs reciprocal lattice vec- 
tor 2M. Thus, the anomalous term of Hsc, which mixes 

and 



XM+k,fc„,tr/; induces a coupling be- 



X-M-k,fc 2 

tween \4>kl) and |<^ kcr ) in H cfi . 

In addition, recall that for a tunneling matrix element 
T between certain electronic states there is a tunnel- 
ing matrix element — T between corresponding charge- 
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conjugated states. This implies that the operator T in 
the Nambu representation Eq. (|5J) is proportional to t z , 
see Eq. 0. 

The single-quasiparticle states are given by the follow- 
ing BdG equations 



(H S c-w)* sc + r*Ti = 0. 



(Al) 



The four-component spinor ^sc (^ti) describes a quasi- 
particle inside the SC (TI). 

Solving the second of the Eqs. (jAlj) for the wave- 
function ^sc and substituting the resultant expression 
in the first one, we derive the effective BdG equation on 
the TI surface 

[H Tl + S(k, w) - ui] *ti = 0, (A2) 
where the self-energy £ on the TI surface reads 

S(k, u) = -j- \T(K k z f r z G sc (k, k„u) t z . (A3) 

2 fez 

The Green's function for the superconductor Gsc = 
(-ffsc — can be derived with the help of Eq. ©: 

G sc (k,fc z ,w) = [e(k,fc z )r 2 + At x - lutq}~\ (A4) 

where e(k, k z ) = (A: 2 + k z ) /2m — e F . Substituting 
Eq. (|A4| in Eq. ([A3]) . one obtains 



/gr A J^A?r,? i^Mi 2 . (A5) 



|A| 2 -uj 2 



Here ro is the unit matrix. The tunneling matrix element 
T(k, k z ) is assumed to vary slowly as a function of k z . 
Then, transforming the k z integral to an energy one, we 
derive 



S(k,o;)«A(k) * T *- UT ° -SUr z 



(A6) 



where A(k) = (ir/2)p(ep , k) |T(k, fc z )| 2 characterizes the 
transparency of the interface. The quantity v{e 1 k) quan- 
tifies the DOS at a given energy and parallel momentum 
in the normal state of the SC. It is equal to 

^±S[e-e(k,k z )}. (A7) 

In what follows we are interested in states close to the 
Dirac cone k = M and ignore the k-dependence of A, 
assuming that A(k) ps A(M). 

The correction to the TI Fermi energy due to tunneling 
is equal to 



6U = 



dk z e\T(k,k z 
2tt e 2 + |A| 2 - 



A(M) / — (A8) 
J-B e e 



A(M)ln|l=0(A), 



where (E e ) is the largest hole (electron) energy rela- 
tive to the Fermi level, and the integral over 1/e is taken 
using the Cauchy principle value. This term is, in gen- 
eral, non-zero. It was discarded in Ref. [l5l presumably, 
because it can be absorbed into the renormalized Fermi 
energy. However, if A depends on the position r (which 
is the case for both our system as well as the system of 
Ref. [l5l ). then 5U becomes spatially inhomogeneous as 
well. Under such circumstances, SU(r) cannot be ab- 
sorbed into the Fermi energy, and has to be treated sep- 
arately. It should be possible, however, to compensate 
8U by an external gate potential. 

Finally, using Hti from Eq. ([5]) and £ from Eq. (jAGD . 
we rewrite Eq. (|A2[) in the explicit forrn^ 



(H eS - u) ti = 0, (A9) 
iv(u)(a-V r )-U(u)) t z + A(w)t x , (A10) 



v^\A\ 2 -uj 2 



y/\A\ 2 - uj 2 + A 



U(oj) 



(U + SU)^^] 2 - uj 2 



A(w) 



^|A| 2 -w 2 + A 
AA 



^|A| 2 -w 2 + A 



(All) 
(A12) 
(A13) 



Appendix B: Bogolyubov-de Gennes equations for 
the radial motion 



In this Appendix we will exploit the cylindrical symme- 
try of the problem, to separate angular and radial vari- 
ables, and to derive the equations for the radial part of 
the wave function. The effective Hamiltonian with a sin- 
gle vortex is equal to 



iv{r) exp (i8a z ) <J x d r -\ — cr y dg 



|A(r)|exp ( i-r z Jr x cxp |~*^ r ' 



(Bl) 
U(r) Tz . 



The first term here is the Weyl-Dirac Hamiltonian in 
polar coordinates. The second term corresponds to the 
anomalous term in the presence of a vortex. 

Equation (|F31|) suggests that it is useful to define a new 
spinor ^ 

*Ti=exp[ifl(r,+<7,)/2]*. (B2) 
Accordingly, the effective Hamiltonian is transformed as 
H' eS = exp [-i6(r z + a z )/2]H cff exp [i6(r z + a z )/2]. (B3) 
The equation for now reads 
1 



a x d r 



-cFyd s + —<J v {t z + a z ) 
r 2r 



IAI 



Ut 7 - u\m = 0. 



(B4) 
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We look for a solution of Eq. (|B4[) in the form \I/ = to ensure single- valuedness of ^ti- It is straightfor- 
gWjj)/^ where \i = 0, ±1,±2, ... is the angular momen- ward to show now that the components of the spinor 
turn. The values of fi are integers (not half-integers) $^ = (/(*, f^, f£, — f£ ) T satisfy Eq. (|22jl . 
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